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ABSTRACT 


We present spectral data cubes of the [Cl] 809 GHz, 115 GHz, 

110 GHz and HI 1.4 GHz line emission from an ~ 1 square degree region along 
the I = 328° (G328) sightline in the Galactic Plane. Emission arises principally 
from gas in three spiral arm crossings along the sight line. The distribution of 
the emission in the CO and [Cl] lines is found to be similar, with the [Cl] slightly 
more extended, and both are enveloped in extensive HI. Spectral line ratios per 
voxel in the data cubes are found to be similar across the entire extent of the 
Galaxy. However, towards the edges of the molecular clouds the [CI]/^^CO and 
^^CO/^^CO line ratios rise by ~ 50%, and the [CI]/HI ratio falls by ~ 10%. We 
attribute this to these sightlines passing predominantly through the surfaces of 
photodissociation regions (PDRs), where the carbon is found mainly as C or C^, 
while the H 2 is mostly molecular, and the proportion of atomic gas also increases. 
We undertake modelling of the PDR emission from low density molecular clouds 
excited by average interstellar radiation fields and cosmic-ray ionization to quan¬ 
tify this comparison. Ending that depletion of sulfur and reduced PAH abundance 
is needed to match line fluxes and ratios. Roughly one-third of the molecular gas 
along the sightline is found to be associated with this surface region, where the 
carbon is largely not to be found in CO. ~ 10% of the atomic hydrogen along 
the sightline is cold gas within PDRs. 


Subject headings: Dark Gas, Molecular Clouds, Carbon, Carbon Monoxide, Galaxy. 
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1. Introduction 


One of the basic activities of a spiral galaxy like our own Milky Way is the continual 
collection of diffuse and fragmented gas and dust clouds into giant clouds of molecules, 


which in turn produce stars (e.g. Molinari et ah 2014). This takes place as part of a cycle 
of matter between the stars and the gas, driven by energy flows arising from radiation and 
mass loss from stars, and the dynamical motions of the gas in the gravitating interstellar 
medium. Cloud formation occurs as the gas cools and the density rises, hrst as atomic 
clouds, then as molecular clouds. Clouds coalesce, cooling is enhanced due to molecule 
formation as radiation is shielded from their interior, and star formation initiated as 
gravitational collapse is triggered inside. 


Determination of the spatial structure and the kinematics of the gas as it transitions 
between the phases of the gas is needed to discern how clouds are formed. While the 
hydrogen gas can be seen from the atomic phase through its principal 21 cm emission line, 
its measurement does not distinguish between the warm and cold phases, only providing 
the total column density of the atomic gas. In the molecular phase the bulk component - 
hydrogen molecules - cannot be seen at all, remaining unexcited in their ground state at 
the typical 10 — 20 K temperatures found in molecular clouds. 

Trace species in the gas are thus needed with emission lines that are sensitive to the 
excitation conditions in the gas, to follow the thermal and chemical transitions that occur. 
After helium (which is inert) and oxygen (whose key lines are blocked by the Earth’s 
atmosphere), carbon is the next most abundant element in the Universe. Carbon can be 
found in ionized (C’*'), atomic (C) and molecular (CO - carbon monoxide) forms in the bulk 
of the dense interstellar medium. All of these species are readily excited in the prevailing 
conditions, with emission lines produced in the terahertz portion of the spectrum for C and 
C"*" (0.5 — 1.9 THz), and the millimeter portion for CO (3 mm = 0.1 THz). 
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CO is readily measured from good observatory sites. C and C’*' are, however, virtually 
unobservable from all but the driest sites on the surface of our planet. We have established 
a new observatory at Ridge A, near the summit of the Antarctic plateau, to open up the 
terahertz spectrum for observation in order to obtain the wide-held, high resolution images 


of the carbon lines needed to pursue this science (Ashley et al. 2012 Kulesa et ah 2013). 


The measurement of the key diagnostic lines from all these species in the atomic and 
molecular phases of interstellar medium we call “following the galactic carbon trail”. 

In this paper we present wide-held data from the Mopra telescope in Australia, 


obtained as part of the Southern Galactic Plane CO survey (Burton et al. 2013 Braiding 


et al. 2015) and the HEAT (High Elevation Antarctic Terahertz) telescope at Ridge A in 
Antarctica of atomic carbon, aimed at pursuing this objective. The data covers roughly 
1 square degree along the G328 sightline through the Galactic plane (i.e. I = 328°). In 
particular, we have examined the data set to determine whether there is evidence for dark 
molecular gas present. This is dehned here as regions exhibiting [Cl] emission but without 
corresponding CO line emission. Practically, since there are no clear regions where [Cl] is 
present but CO is not detected in the moderate beam sizes used (2 arcmin for the [Cl]), we 
search for regions of emission where the [C/CO] abundance may be enhanced. 

Such regions of gas may be expected in the surfaces of molecular clouds. These are 


photodissociation regions (PDRs; e.g. Tielens & Hollenbach 1985 Burton et al. 1990, 1992 


Wolhre et al. 1993), and within extinctions of Ay < 1 mags, from the atomic surface. 


self-shielding may allow signihcant columns of H 2 to exist. However the CO abundance will 
be greatly reduced from the cloud interior, being photodissociated by the far-UV radiation 
that heats and drives the chemistry inside the PDR. Since CO is the normal tracer used to 
indicate the presence of molecular gas, such gas is “dark” to standard survey techniques. In 
these molecular cloud envelopes the carbon will instead be found as either C or C’*', and so 
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be amenable to detection through THz frequency observations. 


Wolfire et ah (2010) model the fraction of such a dark component that exists in giant 


molecular clouds, suggesting that it comprises about one-third of the molecular gas, in 
good agreement with the estimates of the dark gas fraction from gamma-ray observations 


Grenier et al. 

2005 

). In these and previous cloud models (e.g. 

van Dishoeck & Black 

iONN 


the [Cl] emission arises mostly from gas where the hydrogen exists as H 2 but without 
signihcant CO present. Thus [Cl] line emission can be used to trace the dark molecular gas. 


Furthermore, small molecular clouds, here defined as those with total column densities 
such that Ay < Imag. = Nhj = lO^^cm”^, would not be found at all in CO surveys. A 
population of such clouds could remain unseen without a corresponding carbon-line survey 
to detect their presence. 


There is also the possibility of dark atomic gas, where 21 cm HI is optically thick and 
so its line intensity under-estimates the gas column. Analysis of the Planck satellite data 


(Planck Collaboration et al. 2011), comparing the optical depth of the dust emission at 


353 GHz to the column of the HI derived from its 21 cm line emission, suggests that such 
clouds may be widespread across the Galaxy. This can occur when the atomic gas is cold 


(Tg < 80 K), through absorption of background continuum. Fukui et al. (2015) examine this 


further for atomic clouds within a few hundred parsecs of the Sun, but out of the Galactic 
plane. They suggest that, if the column density derived from the dust is indeed all atomic, 
then there are extensive regions where Tg is below 35 K, with optical depths thi as high 
as 3.0; i.e. dark atomic ga^ Garbon would exist as either G or G"*" in such atomic clouds 


^We note that 


Lee et al. 


(2015) argue that these results are incorrect and that the HI is 


not optically thick in the Perseus region that they investigated, though the two studies do 
probe the ISM on very different spatial scales. 
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(but not as CO), and so emission from the THz frequency lines of these species might be 
used as a tracer of such gas. This is especially so in the Galactic plane where there is also 
extensive molecular gas, making it difficult to separate out the dust associated with the 
atomic component. 


By mapping the distribution of carbon along the Galactic plane we can therefore 
estimate where the dark gas may be distributed and how much of it exists within the 
interstellar medium. In this paper we concentrate on an analysis of the carbon emission 
associated with the molecular gas on the Galactic sightline that passes through the G328 
sector, also comparing it the ^^GO, ^^GO, G^®0 and HI line emission. We describe the 
observations in ^and present the data in ^ Examining the variation of [G/GO] at the 
edges of emission features, where the fluxes are smallest, presents several challenges so we 


undertake the analysis through a variety of techniques: flux images (^3.1), histograms of 


line ratio distributions (^3.2), averaged line profiles (|:3.3) and a relatively new technique: 


saturation-hue plots (^3.4). We undertake PDR modelling for the molecular environment 
appropriate for the bulk of the data set in ^ presenting several Figures showing line fluxes, 
optical depths and line ratios as a function of depth into the PDR in §4.1[ We then compare 
these to the data in ® and summarize the results in ® This is Paper II. Paper I (Burton 


et al. 2014) analyzed this same data set, but just for a narrow velocity range associated 


with a single molecular filament. 


2. Observations and Data Reduction 

The data analyzed here comes from three separate spectral line surveys conducted along 
the southern Galactic plane, of THz-band [GI], millimetre-band GO and centimeter-band 
HI line emission. The 62 cm HEAT (High Elevation Antarctic Terahertz) telescope has 
imaged the [GI] 809.3 GHz line (^P 2 —^ Pi) from Ridge A, Antarctica, yielding data cubes 
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with 2' and 0.7km/s resolution (Kulesa et al. 2013). ^^CO, and C^®0 J=l-0 emission 


lines (115.3, 110.2 & 109.8 GHz, respectively) were imaged as part of the 22 m Mopra 
telescope Southern Galactic Plane GO Survey ( [Burton et al. 2013; Braiding et al. 2015), 
with 0.6' and O.lkm/s resolution. The 1.420 GHz HI line comes from archival data in the 


Parkes-ATGA Southern Galactic Plane Survey (the SGPS; McGlure-Grifhths et al. 2005) 
and provides 2' and 3km/s resolution. We will generally refer to these lines as simply [GI], 
^^GO, ^^GO, G^®0 & HI, respectively, in this paper, though sometimes will refer to the [GI] 
line as [GI] 2-1 to distinguish it from the [GI] ^Pi —^ Pq (1-0) line at 492 GHz. We will also 
include discussion of the [GH] ^P 3/2 —^ P 1/2 1.90 THz line, referred to here simply as [GH] 
(and also commonly known as the 158/im line). 

The contiguous data set obtained for the G, GO and H lines shown here covers 


I = 327.8-328.7°, 6 = —0.4 to +0.4°. This data set was first presented in [Burton et al. 


(2014, hereafter Paper I), where an analysis of the emission from a hlamentary feature. 


approximately 75 pc long and 5 pc wide, and with only a narrow 2 km/s FWHM in all 
these spectral lines, was presented, hypothesized to represent a giant molecular cloud in 
the process of formation. In this paper we now examine the data from the complete range 
of velocities covered by these observations with the HEAT telescope, taking in all of the 
emission from across the Galaxy along this sight line. 

The Mopra GO and SGPS HI cubes were re-binned to the same voxel resolution as 
the HEAT GI cubes (both in position and velocity; i.e. 2' and 0.7km/s). All cubes were 
then smoothed with a 2' FWHM Gaussian. The analysis presented here concentrates 
on examining variations seen in the [GI]/^^GO line ratio in the data cubes. The ^^GO 
line is used in preference to the brighter ^^GO line in order to minimize complications 
in interpretation caused by optical depth in the ^^GO line, as the ^^GO line is generally 


optically thin in the Mopra GO survey (see Burton et al. 2013) 


























The full velocity range of the [Cl] emission extends from —120 to 0 km/s. Three 
principal features are seen in its averaged spectral prohle (see Fig. where the [Cl] line is 
shown together with the CO and HI prohles). Fig. shows a schematic model of the Galaxy 


(adapted from Vallee 2014) to aid in the interpretation of these prohles. The three features 
correspond to spiral arm crossings at —100 to —85 km/s, —80 to —65 km/s and —55 to 
—40 km/s, respectively. In turn, they equate predominantly to emission from the Norma 
spiral arm near- and far- crossings, and the Scutum-Crux near-crossing along the G328 
sight line, at distances of ~6, 10 Sz 3.5 kpc from Sun. In addition, the velocity range from 
—80 to —76 km/s corresponds to the quiescent hlament analyzed in Paper I, at ~ 5 kpc 
distance. The data set is split into these four velocity ranges for the analysis described 
below. There is also a weaker feature at —20 km/s arising from the Scutum-Grux far-arm 
crossing, ~ 13 kpc distant. 


A visual rendering of the data cubes for G328 in the four principal lines studied here 
(HI, [GI], i^CO, i^CO) is also shown in Fig. Here the velocity dimension is extended 
along the long axis, and provides a proxy for distance from the Sun (with the proviso of 
the near-far ambiguity). The three features corresponding to the spiral arm crossings are 
readily apparent (Norma near-, Norma far- and Scutum-Grux near-, going from left to 
right). The greater extent of the atomic hydrogen is clear, enveloping the carbon and carbon 
monoxide emitting gas. The carbon also gives the impression of being more extended than 
the carbon monoxide, in-between the distributions for the atomic and the fully molecular 
gas. However care needs to be given towards making such an interpretation, as the relative 
optical depths of the emitting species, as well as the display scales chosen, can result in 
mis-leading visual appearances. In the rest of this Paper we concentrate on a quantitative 
analysis of the relative distributions of the hydrogen, carbon and carbon monoxide emission 
to ascertain whether the impression given here is correct, and to explain why and when 
such a separation can occur. 
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3. Results 

In this section we first consider the integrated flux images for the CO and [Cl] lines 
(i.e. zeroth moment images) in order to compare their morphologies in each of the four 
velocity ranges. While these show suggestions of a more extended distribution for the [Cl] 
emission than the CO this is not conclusive, so we then consider other measures that can 
probe their relative distributions. This includes histograms of the distribution of lines ratios 
and hltering the data to examine averaged prohles from voxels with different values of the 
[CI]/CO line ratio. We also present a new method of presentation (Evans plots) which 
applies a 2D color table (hue and saturation) in order to visualize two variables of interest 
in such a way that pixels with low signal to noise do not dominate a display of the ratio of 
two intensities. 


3.1. Moment Images of Line Flux Distributions 

Integrated flux (Moment 0) images over the four velocity ranges are shown in Fig. 
for both the ^^CO and [Cl] lines, accompanied by an image of the ^^CO overlaid with [Cl] 
contours. These images trace the locations of the molecular clouds within three spiral arm 
crossings along the G328 sight line (Norma far-, Norma near- and Scutum-Crux near-), and 
so are all quite different. The head of the filamentary structure studied in Paper I is to the 
left of centre of the lowest panel in this Figure (~ 328.4°, — 0.1°)[^ 

The ^^CO and [Cl] images clearly look similar in each of the velocity ranges. However 
their overlay in the third column suggests that the [Cl] is, in general, smoother and more 
extended than in the ^^CO images. To check that this conclusion is not an artifact of the 
different original resolutions of the ^^CO and [Cl] data sets the analysis was also repeated 


^The full filament can be seen in CO and HI in Fig. 1 of Paper I. 
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with the cubes smoothed by a 200" Gaussian rather than 120". The same result still holds 
(though is not shown here), of the [Cl] being smoother and more extended than the ^^CO, 
though it is not especially striking when seen at this lower angular resolution. 

This is as expected if the carbon is wrapped around the molecular gas, as well as being 
contained within it. Such behaviour would occur if the outer portion of each molecular cloud 
was lacking in CO. However, while suggestive this result cannot be regarded as definitive 
from these images, so we now examine other metrics to probe the relative distributions of 
the [Cl] and CO emitting gas. 


3.2. Histograms of Line Ratio Distributions 

We show in Fig. [^histograms of selected ratios from the [Cl], ^^CO, ^^CO, C^®0 and 
HI lines in the data. These distributions are determined per voxel (i.e. per 2' spatial pixel, 
0.7 km/s velocity channel). To be included here the data for a [Cl] voxel needed to be > 3 ct, 
whereas for the three CO lines we show histograms for thresholds of la and 3a. These a 
values are determined from the standard deviation of voxels in the continuum portion of the 
data cubes, —120 to —100 km/s. For the HI data, for which a is not well defined since there 
are no clear regions of line-free emission in the spectra, but whose signal-to-noise (S/N) is 
much higher than that of the [Cl] and CO, the threshold was set to an arbitrary value. The 
number of voxels passing these threshold criteria is listed in Table and some statistics 
on the resulting histograms in Fig. [^ are listed in Table [^ These are given for both the 
Icr and 3 ct thresholds applied to the CO fluxes. The mean values and standard deviations 
for the ratios ^^CO/[CI], ^^CO/[CI] and ^^CO/^^CO are 6.8 ± 1.6,1.4 ± 0.5 and 5.5 ± 2.5, 
respectively, when the la threshold for CO is applied. 

The intent here is to search for regions of strong [Cl] and weak (or absent) CO emission. 
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hence the la threshold limit applied to the CO data. There is limited dynamic range in the 
data set, so repeating this analysis with a 3a cnt-off for CO yields few voxels where the 
[CI]/^^CO ratio is “high” (as defined below). We do, however, list metrics in Tablefor the 
corresponding distribntions when a 3 ct threshold for the CO lines is applied. The analysis 
discnssed below yields consistent resnlts with the Icr threshold when this is done, thongh 
from far fewer data points. 


For ^^CO/[CI] the distribution is roughly Gaussian, however for all the other line ratios 
shown in Fig. [^this is not the case. A signihcant non-Gaussian tail exists for larger values 
of the other line ratios. This highlights excess voxels where the line ratio is typically about 
twice the mean value. The tail is especially evident for the [CI]/^^CO ratio, which has a 
mean value of ~ 0.8 and “high” values extending to ~ 2. That the tail is also sensitive 
to the value of the threshold value used to select the CO lines is also apparent from the 
(dotted line) overlays in Fig. where the corresponding distribution for [CI]/^^CO for 
the 3a threshold is almost Gaussian in form. This, of course, means that care is needed 
in interpreting the data in the tail, as noise fluctuations from voxels with weak fluxes will 
cause some ^^CO values to be “low”, and hence the corresponding [CI]/^^CO values to be 
“high”. Nevertheless, by averaging over voxels where the line ratios are high the S/N can be 
improved, allowing us to examine where such voxels tend to be distributed. We undertake 


such an analysis in the following section ((:3.3). 


Fig. § presents a series of scatter plots showing the relationship between the line 
ratios ([CI]/^^CO, [CI]/^^CO and ^^CO/^^CO) and the fluxes of the CO and [Cl] lines, for 
every voxel in the datacube, where the 3a threshold is applied for the [Cl] flux, and la 
for ^^CO and ^^CO. The top panels show the line ratios as a function of the ^^CO and 
^^CO fluxes, the bottom line as a function of [Cl] flux. As these fluxes increase, the line 
ratios tend towards roughly constant values, given by [CI]/^^CO~ 0.5, [CI]/^^CO~ 0.15 
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and i^cO/i^co rNJ 3. In contrast, as the flnxes decrease, the typical values of these line 
ratios tend to increase. Of course, care needs to be given to such an interpretation, for 
noise fluctuations will result in higher derived line ratios for the lowest fluxes when the flux 
threshold applied is also the denominator in the ratio (as it is for the plots in the top line)|^ 
Nevertheless, the median values of the [Cl]/CO lines ratios are clearly rising as the CO line 
fluxes decrease. The overplotted lines in Fig. [^show the weighted minimum least squares 
linear hts to the dataj^ The weak, but clearly negative slopes, to the hts plotted against 
the CO flux, as well as to [CI]/^^CO vs. [Cl], quantify this relationship; e.g. when the ^^CO 
flux is 1 K in a voxel, the [CI]/^^CO ratio is ~ 50% higher than when it is 3K, and for a 
[Cl] flux of 0.5 K in a voxel, the [CI]/^^CO ratio is ~ 40% higher than when it is 1.5 K. 

The ^^CO line flux, being generally optically thin, provides a measure of the product 
Ncofco, the column density of the CO molecules times their Filing factor in the beam, 
whereas for the [Cl] flux it is iVc/cj^ While for any given voxel these factors cannot be 
separated, statistically, as the flux decreases there will be a greater proportion of lower 
column density regions associated with the corresponding voxels. We deduce that for these 
voxels the ratio of the columns of gas in C to CO tends to rise, although the total column 

^Note also that the curvature seen here as the lower boundary in the hrst two plots in 
Fig. [^represents the ratio when the threshold value is applied to the numerator; i.e. the [Cl] 
line. 

^The weighted £t parameters are listed in the caption to Fig. These linear hts are 
strongly weighted towards lower ratios for the lower huxes as the S/N is highest for such 
data points. Hence the £t lines fall below the bulk of the data points. 

^Note that we are implicitly assuming that the [Cl] emission is also optically thin. This 
is consistent with the results of the PDR models we present later, and was also the case for 
the analysis given in Paper I. 
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in both species is lower. 


3.3. Averaged Line Profiles 

We pick a dividing ratio between “high” and “normal” of 1.0 for [CI]/^^CO from 
inspection of the histogram in Fig. and then split the data set into two cubes for those 
voxels that are “high” and those that are “normal”. Visual inspection of the resulting cubes 
suggests that “high” ratio values are generally found around the edges of CO emission 
regions. However, care is needed over any interpretation of high [CI]/^^CO existing at the 
peripheries of molecular clouds because these are also the regions where the S/N is lowest 
since the flux is least. In particular, the la threshold limit applied to the selection of CO 
voxels exacerbates the selection of such voxels at the edges of emitting clouds. 


3.3.1. Splitting the Data based on the [Cl]/^^CO Line Ratio 


We thus apply a mask to separate the voxels with “high” and “normal” ratios, and 
average the resulting cubes to improve the resulting S/N (data from ~ 19, 000 and ~ 55, 000 
voxels are averaged, respectively). The results are displayed in Fig. Here are shown the 
average prohles for ^^CO (in red) and [Cl] (green) in the data cube for all voxels that have 
“high” and “normal” [CI]/^^CO ratios (i.e. > and < 1.0, respectively, as well as passing 


the flux thresholds discussed in [3.2). Prohles for voxels with “normal” ratios are drawn 


with the solid lines, while those for “high” [CI]/^^CO ratios are shown with the dotted 
lines. It is clear that “normal” ratios generally correspond to brighter huxes for both [Cl] 
and ^^CO, as well as having the ^^CO line hux about 50% stronger than the [Cl] hux (i.e. 
red > green). Prohles for “high” [CI]/^^CO voxels, on the other hand, show weaker huxes 
on average, as well as (by design) [Cl] huxes that are (slightly) higher than ^^CO. This 
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result applies across all velocities, and in particular in each of the three principal spiral arm 
crossings along the G328 sightline. Since the resulting prohles have been averaged across 
hundreds of voxels meeting the selection criteria their S/N is good. We can conclude that 
[CI]/^^CO is higher in pixels where the [Cl] and line fluxes are lowest; i.e. as would 

occur preferentially on the edges of molecular clouds. 

There is also the possibility that high ratios occur in voxels which have low optical depth 
(or column density), for, as we discuss in [Q clouds which only extend for ~ 1 — 2mags, 
have a large proportion of their gas in the surface layer (i.e. < 1). Here the carbon will 

predominantly be found as C-atoms (both singly-ionised and neutral). However, as the 
^^CO/^^CO line ratios shown in Fig. [^illustrate, their typical flux ratio is ~ 5 and virtually 
all voxels have [^^CO/^^CO] < 10; i.e. the remains strongly optically thickj^ Thus, 

this possibility (low optical depth) could only apply to relatively few voxels in the data set, 
at most. 


3.3.2. Line Profiles and Ratios as a Function of Velocity 

We now examine the CO, [Cl] and HI lines profiles and ratios as a function of velocity, 
comparing the two groups (“normal” and “high” [CI]/^^CO ratios) with each other. 

Fig. shows the mean prohles for the lines, and Fig. j^for their line ratios, when the 
data has been selected so that the [CI]/^^CO ratio is “high” (i.e. > 1.0) or “normal” (i.e. 

< 1.0) (as well as thresholded, as befor^. The behaviour is similar for the three velocity 
ranges (i.e. spiral arm crossings). The ^^CO and [Cl] huxes are slightly brighter, on average, 

~ 5 when [^^CO/^^CO] = 10, for an isotope ratio [^^C/^^C] of 50. 

^The prohles have been constructed for each voxel meeting the criteria, so the normal- 
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ization is diherent for each velocity channel; hence the diherent appearance from Fig. 
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when the [CI]/^^CO ratio is “normal”, whereas the flnxes are signihcantly stronger for 
“normal” ratios. The HI flux is, however, is slightly smaller for “normal” [CI]/^^CO ratios. 

For the plots showing line ratios (Fig. [^, the [CI]/^^CO ratio is significantly larger for 
the “high” voxels (as it should be, by design). The ^^CO/^^CO ratio is also similarly higher. 
There is, however, little difference in the value of the [CI]/^^CO ratio between the “high” 
and the “normal” voxels. None of these line ratios depend on the velocity, either. They are 
constant, only dependent on the relative value of line flux; i.e. higher [CI]/^^CO and 
^^CO/^^CO ratios for voxels with lower fluxes, and vice-versa. 

The [CI]/HI ratio, however, decreases going from “normal” to “high” [CI]/^^CO values, 
by about ~ 10 — 20%. This indicates that there is slightly more atomic gas associated 
with the “high” [CI]/^^CO ratio gas. Also noticeable is the velocity range for the narrow, 
quiescent hlament (—80 to —76km/s) discussed in Paper I, where the [CI]/HI ratio is 
about 40% higher where [CI]/^^CO is “normal”, indicating substantively more atomic gas 
associated with this hlament than for other molecular gas. 

The conclusion is that high [CI]/^^CO and ^^CO/^^CO ratios occur when the 
hux is weak - i.e. on the edges of the various molecular emission features. At these same 
locations the emission is not, however, signihcantly weaker, on average. This can be 

attributed to optical depth; the line remains optically thick, even when relatively less 

of the species is present, whereas the same sightline samples all of the in the gas, if it 

passes through the edges of clouds. Since the hux is lower here then, as all the CO is 
being sampled, its total column must also be lower. In these same sightlines we hnd that 
the [Cl] hux is only slightly lower, however, and hence so is the column of carbon. Thus, 
elevated carbon to carbon monoxide column density ratios are measured towards the edges 
of all emission features, i.e. at the edges of molecular clouds. 

Furthermore, the lower [Cl]/HI ratios in this same gas indicates that this material is 


16 


associated with relatively more atomic gas at molecular cloud edges. This is as expected for 
sightlines passing preferentially through surface layers {Ay < 1) of PDRs, where CO should 
be photodissociated into C and and a signihcant amount of H 2 is photodissociated into 

H. 

Thus, the relative column of C to CO along these sightlines is increased in comparison 
to sightlines passing through the interiors of molecular clouds. This is also the observational 
signature we have dehned for detecting “dark gas”. This amounts to a statistical detection 
of regions where the carbon (C) is enhanced in abundance compared to the carbon monoxide 
(CO). These must preferentially contain the surface regions of molecular clouds since they 
pass through their edges. Finally, for the quiescent hlament, which we hypothesized to be 
undergoing the process of molecular cloud formation in Paper I, the increased HI abundance 
associated with the cloud edges is consistent with the assertion that it is condensing out of 
the atomic substrate. 


3.4. Hue — Saturation Plots to Examine Variations in Line Ratios 

We present here another technique to examine the signihcance of the increased 
[CI]/^^CO ratios at the edges of the molecular clouds - hue-saturation plots. 

Hue / Saturation plot^ provide a way of visualizing two variables of interest, one of 
which provides some measure of “importance”, and the other its value. This is done by 
assigning the hue of a color to the value and the saturation (or intensity) of the color to the 
“importance”. Hue is cyclic, representing the angle around a color wheel running through 
the spectrum from red to purple (and back to red again); i.e. with range from 0° up to 360°. 

® These are sometimes known as Evans plots; see the description given at 
http://www.ncl.ucar.edu/Applications/evans.shtml, 
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Low importance causes the color to fade into grey regardless of what the value actually 
is. Such plots are useful for examining the ratio of two variables when the overall S/N is 
limited. It avoids the eye being drawn to seemingly high line ratios in regions of low S/N 
data. 


We have constructed such Evans plots to examine the variation of the line ratios by 
using hue (i.e. color) to denote the value of the line ratio, and saturation (i.e. intensity) to 
denote its S/N. Hue is ranged from 0 to 300° in these plots to avoid wrapping, i.e. from red 
for the lowest ratio displayed to purple for the highest. Saturation is ranged from 0 to 1, 
with 1 normalized to the value of the highest S/N. Hue-Saturation-Lightness values (the 
latter, which represents the perceived luminance of the system, is always set to a hxed value 
of 0.5 here) are then converted to red-green-blue for display (using IDL library routines 
which transform between the HSL and RGB color systems). 


Two Evans plots are shown in Figs. 10 - 11 for the velocity ranges of (—80, —76) and 
(—100, —85) km/s; i.e. the hlament and the Norma near- spiral arm crossing (the other two 
spiral arm velocity ranges show similar characteristics and are included in the Appendix; see 


Figs. A.l - A.2). For each velocity range three plots are displayed, [CI]/^^CO, [CI]/^^CO 
and ^^CO/^^CO, so the behaviour of all three lines can be compared. The Evans plots are 
on the left and their corresponding 2D color tables on the right. 


Examination of these images (and in comparison to the flux images of Fig. shows 

that: 


i. For [CI]/^^CO, the highest S/N regions tend to be surrounded by higher-ratio, but 
lower S/N regions. 

ii. For [CI]/^^CO, on the other hand, there are regions of differing line ratio across 
the helds, but with no particular tendency for high ratio regions to surround low 
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ratio regions, or vice-versa. There is also no tendency for lower S/N regions to have 
different ratios than higher S/N regions. 

iii. For ^^CO/^^CO, the line ratios are generally lower within the emission regions where 
there is more CO, and higher at their edges. 

Taking these resnlts in reverse order, we can interpret them as follows: 

iii. Since ^^CO/^^CO variations measnre optical depth changes, we dednce that the 
optical depth is highest (and so the line ratio is lowest) along the sight lines passing 
throngh cores of clonds. It is lowest (with the highest ratios) for sight lines that pass 
throngh the edges of clonds. 

ii. [CI]/^^CO ratio variations trace changes in the relative [Cl] ffnx between molecular 
clouds, since ^^CO is, in general, optically thick and so relatively constant. 

i. Higher line ratios for [CI]/^^CO are seen for sight lines that do not pass through the 
cores of GMCs, but rather tend to pass through their peripheries. Since this ratio 
is sensitive to all the C and CO gas along each sight line, as both lines are largely 
optically thin, then along them the C abundance must be elevated compared to that 
of the CO. This is as expected for the surface layers of molecular clouds, where 
for the first optical depth CO will generally be photodissociated, and so the carbon 
can only be seen as [Cl] (and [CII]), rather than as CO. This is likely mostly dark 
molecular gas, with H 2 existing in the surface layer, and the carbon largely in atomic 
form rather than molecular. 

In the PDR models discussed in the following section the [Cl] peaks in gas which 
is molecular rather than atomic. It may, however, also represent atomic gas, in 
particular if there is [CII] present. As determined in Paper I, the column density of HI 
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is typically ~ 10^^ cm“^ in each voxel. Future measurements, at the same resolution, 
of the 1.9 THz [CII] line, which will co-exist with much of the HI, would allow the 
amount of C’*' and C associated with the atomic gas to be quantified. 


There are also specific regions which stand out in these plots. For instance, at 


(328.65, —0.25) in Fig. 10 both [CI]/^^CO and ^^CO/^^CO are measured to be about twice 
as high as other regions within the hlament, while [CI]/^^CO is unchanged. This can be 
ascribed to a correspondingly lower optical depth in the CO. 


4. Photo Dissociation Region (PDR) Models 


We have examined a range of cloud properties using a modihed form of the Wolfire 


et ah (2010) photodissociation region code to calculate column densities and line emission 


to compare with the data. The model consists of a slab of gas of total optical depth 
Ay,tot illuminated by the interstellar radiation field on two sides. The gas temperature 
and the abundances of atomic/molecular species are calculated as a function of optical 
depth. Ay, under the assumptions of thermal balance, chemical equilibrium and constant 


thermal pressure. For details on the chemistry and thermal processes we refer to Tielens 


& Hollenbach 

li')N5| 

Kaufman et ah 

(2006) 

Wolhre et ah 

(2010 

) and 

Hollenbach et ah 


( 2012 ) 


We have made a number of modifications to the code used in Burton et ah (2014). 


First we have explicitly included chemistry and line transfer in order to model the 
observations. A number of researchers have presented their results on chemistry 


in PDRs including Visser et ah (2009), Rollig & Ossenkopf (2013) and Sziics et ah (2014). 


Here we use the fit by Liszt (2007) to the temperature dependence of the important 


fractionation reaction —)■ -|- 34.8K, and use a ratio of 50 



































20 


as appropriate for the inner Galaxy (Milam et al. 2005). Second, we found that a deeper 


penetration of the external radiation held, compared to our previous models, is required 


to accurately ht the [CI]/^^CO line ratios. Similar results have been reported by Glover 


et al. (2015) and naturally arise in a turbulent gas in which lines of sight pass through 


low column density material. We approximate this process in our models by changing the 
angle dependence of the incident radiation held from isotropic to normally incident. For a 
gas layer at optical depth Ay, the resulting FUV held strength is times larger for 


the normally incident case compared to the isotropic case (see discussion in Wolhre et al 
2010). We used the shielding functions for ^^GO and ^®GO given by Visser et al. (2009) 


which accounts for self-shielding, shielding of ^^GO and ^®GO by H 2 , and shielding of ^®GO 
by i^GO. 

Additional changes involve abundances of metals and Polycyclic Aromatic Hydrocarbons 


(PAHs). Tielens & Hollenbach (1985) noted that an important production reaction for G 


was charge exchange with S (G"*" -|- S —)■ G -|- S+). In addition, it has since been realized 
that ion recombination on PAHs can be important in increasing the abundance of neutrals 
(and also decreasing the abundance of free electrons). We hnd that a decrease in both the 
sulfur abundance and PAH abundance is required to match the ^®GO to [GI] line ratio. 
Our previous models used a sulfur abundance of S/H = 2.8 x 10“^ based on the analysis of 


Savage fc Sembach]| (1996). However, a comprehensive study by Jenkins (2009) hnds a sulfur 
abundance in depleted lines of sight in the solar neighborhood to be S/H = 3.5 x 10“®. 


In a recent study of the [GI] line emission towards the Taurus molecular Gloud, Orr et al 


(2014) hnd a sulfur depletion factor of at least 50 is required to match the observations. In 


light of the more direct observational estimate of the S abundance and the likely presence 


of PAHs, we limit the sulfur depletion to the Jenkins (2009) value. Since we are modelling 


clouds in the inner Galaxy we multiply by a factor of two all gas phase metal abundances, 
dust abundances and PAH abundances to account for the increased metallicity there. This 
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results in a gas phase sulfur abundance of S/H = 7 x 10 ® (see Table [^. 


Even using the depleted sulfur abundance we hud that we still need to reduce the 
Galactic PAH abundance by a factor of two. Reducing the PAH abundance decreases 
the rate of S’*" recombination on PAHs and thus decreases the S abundance in the C’''/C 


zone. This further reduces the neutral carbon abundance. Hollenbach et ah (2012) used 
PAH/H = 2 X 10 ^ for a local abundance. Accounting for higher metallicities in the inner 
Galaxy would result in PAH/H = 4 x 10“^, however we reduce this by a factor of two to 
arrive at a model abundance of PAH/H = 2 x 10“^. There is still considerable uncertainty 
in both PAH rates and PAH abundances and a factor of two decrease is not unreasonable. 

Our “standard” model parameters are listed in Table with some predicted lines fluxes 
and column densities then listed in Table For this standard model we adopt a radiation 
field of Go = 5.1 (with half of this incident on each side; Go is the interstellar held strength 
in units of Habing helds with Go = 1.7 representing an average held strength in the ISM; 


Habing (1968); Draine (1978)). This held is roughly consistent with the interstellar held 


expected at a Galactocentric radius of 5kpc (Wolhre et ah 2003). We also use a constant 


thermal pressure of Pth/k = 2 x 10^ K cm“^. The typical density (when Tgas ~ 30 K) 
is then riHa ~ 600 cm“^. Models are computed as a function of total optical extinction 
through the cloud from Ay,tot = 0.5 mags, for ‘thin’ clouds, up to Ay,tot = 6 mags., where 


Aytot = [A^(H I) + 2A^(H2)]/2.0 x 10^^ cm We demonstrate in (4.1 the ehects of varying 
our standard parameters. 

We have adopted a primary cosmic ray ionization rate of (/crp = 2 x 10“^® s“^, consistent 
with rates for low column density dihuse clouds determined from observations of H^ 


(Indriolo & McGall 2012), OH’*' and HsO"*" (Hollenbach et al. 2012 Indriolo et ah 2015). 


Furthermore, our modelling of the MeV-GeV cosmic-ray huxes (based on Schuppan et 


al. 2012) from the known supernova remnants (G329.7-1-00.4, G327.4-I-01.0, G327.2-00.1, 
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G327.1-01.1, G328.4+00.2 and G327.4+00.4) and Fermi-LAT GeV gamma-ray sources 
(3FGL J1554.4-5315c and 3FGL J1552.8-5330) in the region suggests no enhancement of 
the ionization rate beyond the galactic average (i.e. Ccrp = 2 x 10“^®s“^). There is some 
evidence that the ionization rate is depth dependent in clouds, varying from ~ 2 x 10“^® s“^ 
in low Ay clouds to values of ~ 2 x 10“^^ s“^ in the centers of molecular clouds (e.g. 


As discussed in Paper I, the model has a high cosmic-ray ionization rate but equally good 
hts to the data can be obtained with the low cosmic ray rate. However, we hnd that the 
higher cosmic ray rates produce a slightly higher [GI] line intensity since the higher rates 
allow neutral G to extend deeper into the cloud. 


Padovani et al. 

2009 

Hollenbach et al. 

2012; 

Indriolo & McGall 

2012 

Rimmer et al. |2012 


In addition to providing results from these models for the three principal lines observed 
in the data presented here (i.e. ^^GO & ^^GO 1-0 and [GI] 2-1) we also include the results 
for [GI] 1-0 and [GII] lines emitted at 492 GHz and 1.9 THz, respectively. This is so that 
they may also be applied towards interpreting measurements made with sub-mm telescopes 


in Ghile such as Nanten2, APEX and ALMA (e.g. Graf et ah| 2008 |Lo et ah 2014), as well 
THz frequency telescopes. These latter include airborne ((e.g. SOFIA; [Perez-Beaupuits 


as 


et ah 2012) and space (e.g. Herschel; Velusamy et ah 2010 hanger et ah 2014) facilities as 


well as future telescopes under development for Antarctica such as STO-2 (Walker et ah 


2010 

) and DATE5 ( 

Yang et al. 

2013 


4.1. PDR Model Comparison 


We present here the results for a representative model using our ‘standard’ parameter 
set with Ay^tot = 6 mags. The abundance of G"*", G, ^^GO and ^^GO, as well as of H 2 , as 


a function of the optical extinction. Ay, from the cloud surface is shown in Fig. 12 Both 


the ionized and neutral carbon are seen to be conhned to the surface of the PDR, with the 
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C’*' peaking at the front surface, and mostly disappearing by Ay ~ 0.6 mags, where the C 
peaks. Further in to the cloud the carbon is converted into CO, with rising rapidly 

between Ay ~ 0.5 — 1 mags, followed by a slow rise to Ay ~ 1.8 mags. At this point 
photo-desorption of CO from the cold grains diminishes due to dust attenuation of the FUV 
and freeze-out of the CO onto grain mantles begins, with the CO abundance dropping by a 
factor of ~ 3 at Ay ~ 3 mags, (at the cloud center). The hydrogen gas, on the other hand, 
is converted almost entirely into H 2 by Ay = 0.1 mags, so the C is mostly found in regions 
where the hydrogen is molecular H 2 rather than atomic H. 


The optical depth of the [Cl], [CII] and CO lines as function of Ay in this representative 


model is shown in Fig. 13 The total extinction through the model cloud is Ay^tot = 6 mags., 
and lines here are plotted to cloud center. Thus, the total optical depth in the sightline 
through the cloud is twice that shown at Ay = 3 in the Figure. [CII] rises rapidly until 
Av 0.5 and then remains constant, as by then all the C^ has been converted to C or CO. 
The [Cl] lines rise rapidly to Ay ~ 0.8 past the abundance peak, and then only gradually 
thereafter as the [Cl] abundance falls. The ^^CO continues to rise within the cloud but 
starts to flatten near Ay ~ 2.5 due to ^^CO freeze-out. All the associated lines for these 
three species remain optically thin to cloud center. Through the full cloud extent (i.e. 

Ay = 6) the [CII] and [Cl] 2-1 lines remain thin, while the [Cl] 1-0 and ^^CO lines are 
marginally thick. ^^CO is strongly optically thick beyond Ay ~ 0.8. 


The gas temperature, Tgas, and dust temperature, T^ust, through the FDR are shown 


as a function of the extinction for the representative model in Fig. M The gas temperature 
averages around 30 K within the first magnitude of the extinction, though rises slightly to 
~ 37K for Ay < 0.1, at the very front surface of the cloud. The temperature rise is due 


to molecular hydrogen formation. The collisional rate coefficient of molecular hydrogen 
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exciting is about half that of atomic hydrogen exciting Since the C+ is the 

dominant coolant, the lower excitation rate leads to warmer gas. At deeper layers, the 
temperature falls due to extinction of FUV radiation field resulting in a lower photoelectric 
heating rate. At Ay > 2, the photodesorption of molecules from grain surfaces falls 
sufficiently so that CO stays frozen onto grains. The freeze-out of the dominant coolants 
cause the temperature to rise slightly at cloud center. 


The integrated line intensities for the [Cl], [CII] and CO lines through the PDR as a 


function of column density (i.e. Ay) for the representative model are shown in Fig. 15 The 
intensity plotted is the total intensity emitted through one side of the cloud, including the 
emission from both the near and far sides (i.e., the emission generated at the far side of 
the cloud that passes through the cloud center and emerges at the near side of the cloud). 
These fluxes increase with depth into the PDR, until the carbon has been converted into 
CO, and the CO then frozen-out onto grains. The model line fluxes for the [Cl] 2-1, ^^CO 
and ^^CO lines are typically 6-12 times those measured in the individual voxels. However, 
the model uses a FWHM of 4 km/s, a factor 6 greater than the voxel width of 0.7 km/s 
used in the data analysis, and so is consistent with flux filling factors in the 2 arcmin beam 
of between ~ 50 and 100%. 


In Fig. 16 are shown plots of the [Cl] 2-1 / ^^CO line ratio as a function of optical 


®We note that the collision rates we use for C"*" with H 2 are within a few percent of the 


recent rates of para-H 2 with C"*" in Wiesenfeld & Goldsmith (2014). Adopting these rates 
results in only a few percent decrease in the gas temperature. We also note that our collision 


rate of H with C"*" (from Flower 1990) is slightly higher than that in Barinovs et ah (2005). 


Adopting this more recent rate would further increase the gas temperature in the atomic 
H layer but since collisions with atomic hydrogen dominate in only a narrow region at the 
outer edge of the cloud this makes a negligible difference to our results. 
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extinction, Ay tot, through molecular clouds for a variety of values of the input parameters. 
For each curve we keep the parameters fixed for the ‘standard’ model, and vary only the 
parameter that is listed in the figure legend. The PAH/H curve is for a PAH abundance 
that is twice as high as for the standard model while the S/H curve is for a sulfur abundance 


that is twice the abundance used in our previous modelling (Kaufman et ah 1999) and 


eight times larger than in our standard case. The Go curve is for a radiation held that 
is the local Galactic held, and the thermal pressure curve is a factor of 3 lower than the 
standard model. The parameters range roughly between those expected for nearby clouds 
and those in the inner Galaxy. The higher PAH/H and S/H curves increase the [GI]/^^GO 
ratio by increasing the G abundance as described in ^ The lower Go model produces a 
lower [GI]/^^GO ratio due to a lower heating rate in the [GI] region. For the low pressure 
models, the lower density pushes the molecular transition deeper into the cloud, thereby 
reducing the ^^GO column for the low Ay^tot models. 

The line ratio is seen to decrease as the cloud column density increases, since the [GI] 
emission is confined to the front surface. Ay < 0.6 mags, (see Fig. 15), whereas the ^^GO 


emission extends deeper into the cloud. Except for the low thermal pressure model, beyond 
Aytot ~ 4 mags, the line ratio remains constant, indicating that neither line is contributing 
signihcant additional line flux for thicker clouds. All models obtain the observed [GI]/^^GO 
ratios of < 2 for sufficiently thick clouds. The lowest allowed column density cloud is 
Aytot ~ 1, for the low FUV field model. 

The dependence of the ^^GO/^^GO line ratio as a function of the total optical depth 


through the cloud, Aytot is shown in Fig. 17 Again, the line ratio is high for low optical 


depths, but rapidly decreases as Ay rises, plateauing beyond ~ 4 mags., similarly to 
[GI]/^^GO. In this case, however, this is caused by the freeze-out of the GO molecules onto 
dust grains. 
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The dependence of these two line ratios (^^CO/^^CO and [Cl] 2-1/^^CO) are also 
shown as a function of the line flux in Figs. 18 and 19 Similar behaviour to the plots 


against extinction, Ay, are exhibited, with the line ratios being high for small values of the 
line flux. This reflects, correspondingly, the CO optical depth being lower for thin 
clouds (i.e. Ay < 2), and the [Cl] emission arising from the same region, near to the cloud 
surface. In practice, however, observing such high ratio but low intensity gas in our data 
set will be difficult since one or both lines may not be detectable. We also caution that low 
intensity gas could have low ratios, if the beam Filing factor is small. In principle, however, 
it is possible to have gas with ^^CO/^^CO and [CI]/^^CO ratios much higher than those 
observed in the G328 data set presented here. They would also be more readily detectable 
with higher angular resolution measurements, able to separate such regions out from high 
column density sightlines. 

We also calculated the ‘dark gas’ contribution in the standard model. Using the 
dehnitions in Wolhre et al.| ( |2010 ), we assume the dark gas layer begins when the gas is 
half molecular (i.e., 2n(H2) = 0.5n) and ends when the optical depth in the CO (1 — 0) 
transition equals one. We hnd the line-of-sight dark gas fraction is 27% of the total column 
density of this = 6 mags. slab. Since the H 2 forms essentially at the cloud surface 

(dlv 0.1), the dark gas fraction of the total molecular column density is also about 27%. 
We can also estimate the fraction of dark gas in a spherical cloud by converting the Ay 
steps to radius steps (Sr = SAy(r) 2.0 x 10^^/n(r), where Sr and SAy are the distances 
between radius and Ay grid points, respectively). The mass at radius r = T,Sr is found by 
summing the density in spherical shells M(r) = Sn(r)47rr^5r (H cm“^). We hnd a fraction 
of dark gas compared to the total molecular gas is ~ 67% with a similar fraction when the 
total (molecular plus atomic) mass is included. The dark gas fraction is a factor of two 


higher than that found by Wolhre et ah (2010) for local Galactic GMGs. However, the 


mean area-averaged column density for the model cloud is found to be Ay = 2.7. This is 
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smaller than that for GMCs [Ay ~ 8; Solomon et ah| 1987) and the dark gas fraction is 


predicted to increase as the mean column density decreases (Wolhre et ah 2010) 


Finally, we have estimated the optical depth of the (cold) atomic gas in the PDR 
through application of the radiative transfer equation N = TsAVX{1 — e“^), where N is 
the column density of atomic hydrogen in the PDR model (Table |^, Ts its spin temperature 
(i.e. as at the front of the PDR; 35K from Fig. 14), AV the FWHM line width (Table |^, 


and X the HI X-factor (1.8 x 10^® cm ^ K ^ km ^ s; Dickey & Lockman (1990)). This yields 


r ~ 0.3 for the atomic gas in the PDR. Following Goldsmith & hanger (1999), the antenna 


temperature for the 21cm HI line emission from this gas is then given N = Ta AVX , 
which yields Ta ~ 8 K. This can be compared to the typical brightness temperatures 
observed along the G328 sightline, 90-100 K (Fig. [^, which are dominated by hydrogen 
in (warmer) atomic gas clouds rather than in PDRs. From their ratio we estimate that 
the fraction of atomic hydrogen found in PDRs along the sight line is ~ 10%, with the 
remaining ~ 90% being in atomic clouds. Furthermore, we may also estimate the amount of 
“dark” atomic gas in the PDR by calculating the column derived by applying the optically 
thin limit in the above formula, and comparing this to the actual column from the model. 
We hnd that ~ 15% of the atomic gas in the PDR is dark; i.e. a column of H amounting to 
~ 10^®cm“^. This is also a negligible fraction of the total column of atomic gas along the 
G328 sightline, around ~ 1% of it. 


5. Discussion 


The models presented above in 14.1 show that, in the low-FUV field PDR environment 
typical of the diffuse molecular environment, that ionized G"*" is found within Ay < 0.5 mags, 
of the cloud surface. Neutral G extends as far as Ay ~ 2 mags, into the PDR, but peaks 
in abundance at Ay ~ 0.5 mags. Molecular GO exists deeper in than Ay ~ 0.5 mags. 


























Freeze-out of CO begins by Ay ~ 2 mags, with only one-third of the molecule remaining in 
the gas-phase by Ay ~ 3 mags. 

Given the different extents where these species are found inside a PDR, the measured 
line ratios, such as ^^CO/^^CO and [CI]/^^CO, will depend on the column density of the 
PDR sampled by the line of sight. In particular, those sight lines that either only pass 
through their edges, or include low column density PDRs (i.e. only extending for 1-2 mags, 
rather than the more typical 8 mags, for GMCs) will display elevated values of these ratios. 

Optical depth, however, also needs to be accounted for in determining line fluxes, and 
hence line ratios. [GI], [GII] and ^^GO remain optically thin within these emitting regions. 
However, the ^^GO line rapidly becomes optically thick further into the molecular gas from 
the position where it is first encountered. 

Thus [GII] line flux rise until a depth of Ay ~ 0.5 mags, is reached into a PDR (and 
then remains flat as the sightline penetrates further), while [GI] rises until a depth of 
Ay ~ 1 mags. In both cases a flattening in flux with further depth is because the species 
are not found further into cloud. For the ^^GO line flux, in contrast, while it rises rapidly 
between Ay = 0.5 and 1 mags., it then saturates as it becomes optically thick. Only modest 
rises in flux then occur going deeper. On the other hand, ^^GO line fluxes continue to rise 
steadily with increasing depth into a PDR until GO freezes out. Representative line fluxes 
and column densities, for a PDR with a total column of Ay = 6 mags., are listed for the 
standard model parameters in Table The column of G and G"*" are both approximately 
one-third that of the total column of gas-phase GO found though the PDR. 

Gomparison with the flux measurements in Fig. shows that these model predictions 
are consistent with the data for ^^GO, ^^GO and [GI] 2-1, if the average beam hlling factor 
in the 2' voxel is ~ 50% (and taking into account the difference between the model 4 km/s 
FWHM and the voxel channel width of 0.7 km/s). The corresponding [GI] 1-0 and [GII] 
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fluxes are thus predicted to be of order 7 and IKkm/s, respectively, if observed with the 
same spectral and spatial resolution. 

The models show that the [Cl] 2-1/^^CO and ^^CO/^^CO line ratios are constant for 
clouds with total extinctions Ay,tot >2 — 4 mags, (dependent on the particular model), with 
ratios of ~ 1 and ~ 4, respectively. For lower Ay,tot these line ratios rise, as the relative line 
fluxes are sensitive to the depth of the PDR. Small Ay,tot corresponds to PDRs which are 
either relatively thin, or for sight lines which are predominantly passing through the surface 
layers of PDRs. This latter case would also correspond to sight lines passing through the 
edges of molecular clouds rather than through their cores. The ratios can rise to high values 
(> 3 and > 15, respectively) for small Ay,tot, but then the corresponding line fluxes are 
themselves small as only a limited column of gas is being sampled. No data with such high 
ratios was seen, though possibly with higher angular resolution (e.g. as with ALMA), when 
beam dilution of any such smaller clouds that might exist would be much reduced, then 
higher ratios might be found in some gas. 


On the other hand, limb brightening would result in larger line fluxes. Since the 
measured line fluxes are actually smaller when [CI]/^^CO is larger (see Figs. |^and[^, 
this implies that any limb brightening that occurs is overwhelmed by smaller beam filling 
factors. 


The typical [CI]/^^CO line ratios measured here are ~ 0.7, rising to rsj . 3 (Fig.i. The 
latter are found where the fluxes are, in general, weaker (Fig. |^. While for individual pixels 
the effects of weak intrinsic flux or beam dilution cannot be disentangled, the modelling 


shows that if the is intrinsically weak then the [CI]/^^CO ratio is higher (Fig. 19). 
This is consistent with the associated data being dominated by voxels with intrinsically 
weaker flux rather than by beam dilution. 

5, rising up to 15 (Fig.g. As for 


Typical ^^CO/^^CO line ratios observed are 
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[CI]/^^CO, the latter generally occur where fluxes are lower. And similarly (Fig. 18), 
these values are consistent with the higher ratios being associated with sight lines passing 
through lower column density gas (and/or just the front surfaces of PDRs). 

To summarize, the line fluxes for ^^CO, and [Cl] presented here for the G328 

sightline of the Galaxy are consistent with being produced in low density (riHa ~ 600 cm“^) 
photodissociation regions being excited by far-UV radiation fields with strengths around 
the typical average interstellar value (Gq ^ 1.7 — 5.1). Larger line ratios of [CI]/^^CO are 
found, in general, associated with smaller values of the line flux. These smaller line 

fluxes are, in turn, generally associated with being emitted by gas with lower intrinsic 
line fluxes, rather than arising in sight lines with more beam dilution (which is found to 
average around 50%). Such regions are associated with the front surfaces of PDRs, within 
an optical extinction of 2 mags, from the atomic hydrogen interface. Here carbon can be 
found in ionized (C’*'), neutral (C) or molecular (CO) form, rather than being dominated by 
CO at greater depths. Such regions are also associated with dark molecular gas, molecular 
clouds (i.e. H 2 ) where CO does not provide the dominant tracer for the gas. 


6. Summary 

We have presented data cubes showing the distribution of the [Cl] 2-1, ^^CO, ^^CO 
and C^®0 1-0, and HI line emission, over a ~ 1° region of the Galactic plane along a 
sightline towards G328, with angular and spectral resolutions of ~ 1' and ~ 1 km/s. The 
[Cl] data comes from a new telescope, HEAT, sited on the summit of the Antarctic plateau 
at the driest location on the Earth, where the THz atmospheric windows are opened for 
observation. The CO data was taken with the Mopra telescope, and the HI data with the 
Parkes and ATCA telescopes, all located in Australia. 
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Complex morphology is evident in all species, with the [Cl] and CO line emission 
extending over 120 km/s in extent, and arising principally from molecular cloud complexes 
in three spiral arm crossings along the sight line. The distribution in these atomic 
and molecular species is very similar in both angular and spectral dimensions, and is 
encompassed by more extensive HI line emission. However close examination of the [Cl] 
and CO emission shows the latter to be slightly more extended spatially. 

Average line ratios for ^^CO/[CI], ^^CO/[CI] and ^^CO/^^CO are found to 
6.8 ± 1.6,1.4 ± 0.5 and 5.5 ± 2.5, respectively, across all voxels in the data cubes when 
resampled to the same resolution. These ratios are relatively constant with velocity, across 
the three spiral arms crossed along the sight line. However, at the edges of the emission 
features the [CI]/^^CO and ^^CO/^^CO ratios are typically found to be 50% larger, with 
[Cl]/HI ratios around 10% lower. This is attributed to relatively more C than CO at 
the edges of molecular clouds, with the ^^CO intensity providing the best measure for 
the column of CO, rather than the optically thick ^^CO. On the other hand, the relative 
amount of atomic gas is seen to rise at the edges of the clouds. 

PDR models were constructed for diffuse molecular gas exposed to average interstellar 
radiation fields to explore the behaviour of the ^^CO, ^^CO, [Cl] 1-0 and 2-1, and [CH] line 
emission with increasing extinction. Ay, into molecular clouds from their front (atomic) 
surface. Charge exchange reactions between C"*" and both S and PAHs were also found to 
be important to match the data, and require signihcant depletions of the S and PAHs. The 
models then reproduce the broad behaviour seen in the data, in particular the line fluxes 
(with average beam hlling factors of ~ 50%) and the increased line ratios for [CI]/^^CO 
and ^^CO/^^CO in the surface layers, i.e. when Ay < 2 mags. The highest values of the 
[CI]/^^CO ratio measured require sightlines passing through clouds that have extinctions 
Ay < 1 from their exciting sources. These may arise from sightlines that only pass through 
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the edges of clouds, and/or small clouds whose total column density extends for Ay ~ 1. 
The and [Cl] lines remain marginally optically thin through the PDR, and so provide 
direct probes of the amount of CO and C present, respectively. 

Roughly one-third of the molecular gas along the sightline is estimated by these models 
to be associated with dark molecular gas, and about two-thirds in total when accounting 
for a spherical geometry. This is gas where CO line emission does not provide a full tracer 
for the column of molecular gas that is present. Considering the column of cold atomic gas 
in the front surface of the PDR, we also estimate that it contains ~ 10% of the total atomic 
gas along each sight line (the rest being in atomic gas clouds). Of the atomic gas within 
the PDR, ~ 15% of it is dark, i.e. the fraction not measured through use of the standard 
optically thin assumption for determining column densities. This is only ~ 1% of the total 
atomic column along sight line, however. 

The modelling would also be aided by measurements of the [Cl] 1-0 line emission at 
492 GHz, whose ratio with the [Cl] 2-1 line we have observed at 809 GHz provides a probe 
that is sensitive to the gas temperature, predicted here to be ~ 30K for Ay < 1. Such 
measurements could be made from the Atacama plateau in Ghile, where the the superior 
angular resolution obtainable with larger telescopes such as Nanten2, APEX and ALMA 
would allow for a better comparison of the relative distributions of the G and GO, and so 
in determining more clearly the form and extent of the carbon envelope around molecular 
clouds. 


A full quantihcation of the amount of dark molecular present also requires measurement 
of the [GH] line emission, which is conhned within Ay < 0.5 of the PDR surface. Some 


airborne and space measurements currently for such [GH] now exist (e.g. Pineda et ah 


2014 Velusamy et ah 2015). We provide predictions of its intensity and distribution with 


extinction into molecular clouds, to aid in their interpretation and for future observations 
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made with THz facilities in Antarctica, such as the STO-2 long duration balloon-borne 
telescope scheduled for launch in 2016 and China’s proposed 5m DATE5 telescope for Dome 
A on the summit of the Antarctic plateau. 
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Table 1. Statistics on the Voxels used for the Histogram Analysis 


Line 

Number 

na 

Threshold (K) 

Number 

na 

Threshold (K) 


204,350 

1 

0.9 

90,059 

3 

3 


148,623 

1 

0.3 

38,379 

3 

0.8 

Q 

00 

O 

91,875 

1 

0.2 

2,783 

3 

0.6 

[CI] 

83,857 

3 

0.4 




HI 

457,746 

Arbitrary 

40 





Note. — Number of voxels (i.e. pixel-velocity channels) that pass the threshold criteria 


used for the histogram ratio analysis (see (:3.2), together with the number of standard devi¬ 
ations (cr) these correspond to and the corresponding threshold values (in K). For the three 
CO lines, the corresponding number of voxels that pass the threshold criteria if 3 ct is used 
instead of la is also listed. The HI limit is arbitrary since there are no line-free channels to 
estimate the standard deviation from. However, the S/N in the HI is very much higher than 
for the other lines, so this limit has little effect on the results. 
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Table 2. Line Ratio Statistics 


Ratio 

Pair 

# 

Mean Med. Max. 

Icr threshold for CO 

a 

# 

Mean Med. Max. 

3ct threshold for CO 

a 

12C0/[CI] 

83,545 

6.8 

6.7 

18 

1.6 

72,112 

7.1 

7.0 

18 

1.5 

[Clj/i^CO 

83,545 

0.16 

0.15 

0.7 

0.04 

72,112 

0.15 

0.14 

0.34 

0.03 

i^CO/i^CO 

115,687 

5.5 

5.0 

25 

2.5 

36,301 

4.3 

4.2 

10 

1.2 

[Clj/i^CO 

74,063 

0.84 

0.76 

3.0 

0.33 

36,751 

0.63 

0.61 

1.6 

0.14 

13C0/[CI] 

74,063 

1.4 

1.3 

4.1 

0.5 

36,751 

1.7 

1.6 

4.1 

0.4 

12C0/C1®0 

46,750 

10 

8.7 

62 

6.1 

1,715 

7.8 

7.0 

22 

3.2 

12CO/HI 

203,197 

0.03 

0.02 

0.34 

0.02 

89,828 

0.05 

0.04 

0.34 

0.02 

[Cl]/HI 

83,588 

0.007 

0.006 

0.047 

0.003 






i^CO/HI 

142,961 

0.007 

0.005 

0.12 

0.006 

38,181 

0.007 

0.012 

0.12 

0.008 


Note. — Statistics for the voxels which exceed the threshold limits in Table [T] for each 
line pair. Columns 2-6 refer to the la threshold for the CO lines, columns 7-11 to the 3 ct 
threshold. is the number of voxels that meet the threshold limits. Both the mean and 
median values are shown, as well as the maximum value measured for the ratio and the 
standard distribution of the distribution. For [Cl]/HI only one set of threshold limits has 
been used, hence the values are only tabulated once. Fig. shows the histograms for the 
corresponding line ratio distributions. 
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Table 3. PDR Model Inputs 


Input Parameters for the Standard Model 

Parameter 

Value 

Go 

5.1 Habings 

Ay,tot 

6 mags. 


(unless otherwise specihed) 

Pttik 

2.0 X 10^ K cm-3 

AR (FWHM) 

4.0 km s“^ 

Ccrp 

2.0 X 10-1® s-i 

[12Q/13Q] isotopologue ratio 

50 

[C/H] abundance 

3.2 X 10-^ 

[S/H] abundance 

7.0 X 10-® 

[PAH/H] abundance 

2.0 X 10-^ 


Note. — Parameters for the 2-sided PDR models (i.e. the model is ID with radiation 
incident from both sides). Gq is the radiation field strength in free space in units of the 
Habing held (1.6 x 10“^ ergss“^ cm“^; a representative value for the average interstellar 
radiation held). Pth/k is the thermal pressure in Kcm“^. AV is the FWHM line width 
in km s“^. Ccrp is the primary cosmic ray ionization rate per H nucleus. 
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Table 4. PDR Model Outputs 


Output Predictions from the Standard Model 
Parameter Value 


IC^CO 1-0) 

50 K km/s 

IC^CO 1-0) 

12 K km/s 

/([Cl] 1-0) 

14 K km/s 

/([Cl] 2-1) 

5 K km/s 

/([CII]) 

2 K km/s 

/(HI) 

32 K km/s 

N(i2C0) 

9.5 X lO^'^cm”^ 

N(i3CO) 

2.1 X lOi^cm-2 

N(C) 

3.3 X lOi^cm-2 

N(C+) 

2.9 X lOi^cm-2 

N(H) 

7.1 X 10^® cm“^ 

Tgns 

18 K 

Tdust 

9 K 

^gas('rcO = 1) 

28 K 

/dust ('Too = 1) 

14 K 

^H2('rcO = 1) 

580 cm“^ 
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Note. — Output parameters from the standard model of Table are for a PDR 
extending to an optical extinction Ay,tot = 6 mags. (i.e. Ay = 3 to its center). Column 


densities (in cm ^), temperatures (in K) and line fluxes (in K km/s) are as per Figs. 12 


14 and 15, with the latter including the contributions from the near and far sides of the 


cloud. Also given are the temperatures and H 2 density when tqo = 1- The fluxes are for 
the model FWHM of 4 km/s, and so need to be divided by 6 to be directly compared 
with the data, for which the voxel width is 0.7 km/s. 





Figures 
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Fig. 1.— Line profiles of the mean integrated emission from the entire region covered by 
the ~ 1° X 1° region of the HEAT G328 data cube. Lines are ^^CO/5 (blue), (red), 

(magenta), [Cl] (green) and HI/150 (black). Three principal spiral arm crossings can 
be seen in the CO and [Cl] data between —100 to —85 km/s (Norma; near-portion of arm), 
—80 to —65 km/s (Norma; far-portion of arm) and —55 to —40 km/s (Scutum-Crux; near¬ 
portion), respectively. The weaker feature at —20 km/s corresponds to the far-portion of 
the Scutum-Crux arm. 
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Fig. 2.— Schematic of a four-arm model of the Galaxy to aid in visualization of the data 
along the G328 sight line. The model uses the parameters from [Vaflie (2014) with the spiral 
arms shown by color lines: Perseus (turquoise), Sagittarius (green), Scutum-Grux (red) and 
Norma (black). Their extent is indicated by the corresponding dashed lines. The G328 
sightline from the Sun is indicated by the wedge, with the Solar Gircle (at Rq = 8kpc) 
shown by the dashed purple line (the cross indicates the Galactic center). The color shading 
within the wedge shows the expected line of sight radial velocities (14sr) using the galactic 
rotation model of McGlure-Grifhths & Dickey (2007[) for the Fourth Quadrant, matched to 


the Brand & Blitz (1993) model for the outer Galaxy (i.e. positive velocities, when R > Rq^ . 


The short orange arc indicates the locus for the tangent point (where radial velocities are at 
their most negative). The spatial scale along the axes is in kpc, and the numbers along the 
wedge show the distance from the Sun, also in kpc. 
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12CO 13CO 



Fig. 3.— Renderings of the G328 data cnbes for (bine), (red), [Cl] (green) and 
HI (yellow). Galactic longitude and latitude are along the short axes, with velocity along 
the long axis (with the most negative velocities to left). The three principal arm crossings 
along the sight line can be distinguished: Norma near-, Norma far- and Scutum-Crux near-, 
going from left to right. These are also shown as moment images in Fig. The relative 
extent of the atomic and molecular gas can be gauged, with the atomic hydrogen enveloping 
both the carbon and carbon monoxide emitting-gas. 




Galactic Longitude Goloclic Longitude Galactic Longitude 






Galactic Longitude Galactic Longitude 


G3 S6 v=—60 to —76 km/s. 13C0 image, [Cl] contou rs 



336.B” 326.4' 326 2' 326 0' 327 8' 

Galactic Longitude 


Fig. 4.— Across: (i) Column 1: moment 0 images (i.e. integrated flux), overlaid with 

their own contours (in black, at every 10% of the peak value), (ii) Column 2: [Cl] moment 
0 images, overlaid with their own contours (in red, at every 10% of the peak value) and 
(hi) Column 3: ^^CO images overlaid with red [Cl] contours. Down: Velocity ranges for (a) 
—100 to —85 km/s (Norma near), (b) —80 to —65 km/s (Norma far), (c) —55 to —40 km/s 
(Scutum-Crux near) and (d) —80 to —76 km/s (the filament studied in Paper I). The scale 
in the color bar refers to the flux in image in K km/s. The original datacubes were all first 
smoothed with a 120" FWHM Gaussian. 
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Fig. 5.— Normalized histograms of the distributions of selected line ratios from the data set 
(^2CO/[CI], [CI]/^2co, i2co/^3co, i3co/[CI], [CI]/i3CO, ^2C0/C^^0, ^2C0/HI, [CI]/HI& 
[CI]/^^CO, respectively, as labelled) determined on a voxel-by-voxel basis (i.e. pixel-velocity 
channel). To be included the [Cl] flux needed to be > 3 ct per voxel but the fluxes for each 
of the 3 CO lines (^^CO, ^^CO, C^®0) only needed to be > Icr in the voxel (see Table [^. 
Overlaid as dotted lines are the corresponding distributions if a 3a threshold is applied for 
the three CO lines instead. Statistics relating to the distributions are listed in Table 
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Fig. 6.— Distributions of various line ratios vs. line intensities, la thresholds are applied 
to CO, 3a for [Cl] (see text). Shown are, on top row, then second row: (a) [CI]/^^CO 
vs. ^3co, (b) [CI]/i2CO vs. 12CO, (c) i^co/i^CO vs. ^^CO, (d) [Clj/^^CO vs. [Cl], 
(e) [CI]/^^CO vs. [Cl] and (f) ^^CO/^^CO vs. [Cl]. Ten adjacent points have been av¬ 
eraged together in making these plots in order to improve their clarity and error bars 
have also not been shown. The lines show the best linear fit (weighted, minimum least 
squares) to the data; the coefficients (A, B), with Ratio = A x Flux -|- B are as follows: 
(0.8,-0.1), (0.2,-0.004), (5,-0.7), (0.7,-0.2), (0.1, 0.02) and (5,-1.4), respectively [and, as 
discussed in the text, error bars are lower for the smallest ratios when the flux is low, hence 
the linear fits are weighted to pass through these points rather than the bulk of the points, 
which have lower S/N]. The number of points fitted in each plot is evident from Table 
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Fig. 7.— Averaged line profiles for and [Cl] over the entire data cube for all voxels 

that have “high” (> 1.0) and “normal” (< 1.0) [CI]/^^CO line ratios. is shown in 

red and [Cl] in green. Voxels with “normal” ratio values are the solid lines and those with 
“high” ratios are dotted. Clearly those regions with “high” ratios also have lower fluxes, on 


average. 
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Velocity Vlsr (km/s) 

Fig. 8.— Mean line brightness per velocity channel for the ^^CO/5 (bine), (red), [Cl] 

(green) and HI/150 (black) lines. Solid curves are for voxels where [CI]/^^CO <1.0 (i.e. 
“normal”) and dotted for voxels where [CI]/^^CO > 1.0 (i.e. “high”). The velocity ranges 
shown are chosen for the three spiral arm crossings discussed in (i.e. —100 to —85 km/s, 
—80 to —65 km/s and —55 to —40 km/s, with the narrow hlament being contained within 
the hrst of these). The flux in a voxel needs to exceed the thresholds in Tablej^to be included 
in this analysis. Prohles are then averaged over all the voxels which exceed the thresholds, 
and placed into the appropriate line ratio range (i.e. “high”, “normal”) for display. Voxels 
that do not meet the threshold criteria are excluded from the averages shown here. The 
normalization for each velocity channel is thus different than in Figs. [T]and[^ 
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Velocity Vlsr (km/s) 

Fig. 9.— As for Fig. but showing the averaged line ratios for (i)[CI]/^^CO (in red), 
(ii)[CI]/^^CO (in blue), (hi) [^^CO/^^CO] /lO (in magenta) and (iv) [CI]/HI xlOO (in black), 
instead of the line intensities, for the same velocity ranges. As before, solid curves are for 
voxels where [CI]/^^CO < 1.0 (i.e. “normal”) and dotted for voxels where [CI]/^^CO > 1.0 
(i.e. “high”). It is clear that for the “high” voxels the [CI]/^^CO and ^^CO/^^CO ratios 
nearly double, [CI]/^^CO is unchanged and [Cl]/HI decreases by ~ 10% (and noticeably 
more so in the velocity range of the narrow hlament, —80 to —76 km/s). 
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Fig. 10.— Evans plots for the velocity range —80 to —76 km/s corresponding to the hlament 
(left) and their corresponding 2D color tables (right). Shown are, from top to bottom, 
[CI]/^^CO, [CI]/^^CO and ^^CO/^^CO. Color (line) denotes the valne of the line ratio x—axis 
of the color table), and the satnration (intensity) its S/N, with low S/N pixels fading to grey 
(y—axis of the color table). 
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Fig. 11.— Evans plots for the velocity range —100 to —85 kni/s corresponding to the Norma 
near-spiral arm crossing (left) and their corresponding 2D color tables (right). Other details 


are as for Fig. 10 
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Fig. 12.— Abundance of H2/10^ (black dash), C’*' (yellow dash-dot-dot), C (green dot- 
dash), (blue long-dash), and ^^COxlO (red solid), as a function of optical depth Ay 

(in magnitudes) for the standard cloud model. 
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Fig. 13.— Optical depth of the emission lines for [Cl] 1-0 (brown) and 2-1 (green), ^^CO/15 
(blue) and (red) 1-0, and [CII] (yellow), as a function of Ay for the standard model. 

The total extinction through the cloud is Ay^tot = 6 mags., and lines are plotted to the cloud 
center. The total optical depth in the line through the cloud is thus twice that shown at 
Ay = 3. 
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Av 

Fig. 14.— The gas temperature, Tgas, and dust temperature Tdust as a function of Ay for 


the standard model. 
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Fig. 15.— Integrated intensities, in K km/s, of the emission lines for ^^CO/3 (blue), 

(red) 1-0, [Cl] 1-0 (brown), [Cl] 2-1 (green) and [CII] (yellow), as a function of optical 
extinction. Ay, into the PDR, for the standard model. The total extinction in magnitudes 
through the cloud is Rytot = 6. The intensity plotted is the total intensity emitted through 
one side of the cloud including the emission from both the near and far sides. 
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Fig. 16.— The [Cl] 2-1 / 1-0 line ratio as a function of total optical extinction through 

the cloud, Ay tot, for a range of model parameter values, as indicated in the legend. 
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Fig. 17.— The ^^CO/^^CO J=l-0 line ratio as a function of total optical extinction through 
the cloud, Ay^tot) for a variety of model parameters as indicated in the legend. 












co(i-o)/’^co(i-o) 


62 



Fig. 18.— ^^CO/^^CO J=l-0 line ratio as a function of the line flux for a variety of 


model parameters, as indicated in the legend. 
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Fig. 19.— [Cl] 2-1 / line ratio as a function of the line flux for a variety of model 


parameters, as indicated in the legend. 











64 


A. Appendix 


Two Evans plots are shown in Figs. A.l - A.2 for the velocity ranges of (—80, —65) and 


(—55, —40) km/s; associated with Norma far- and Scutnm-Crux near- spiral arm crossings. 
For each velocity range three plots are displayed, [CI]/^^CO, [CI]/^^CO and ^^CO/^^CO, so 
the behavionr of all three lines can be compared. The Evans plots are on the left and their 
corresponding 2D color tables on the right. 
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Fig. A.l.— Evans plots for the velocity range —80 to —65 km/s (left) and their corresponding 
2D color tables (right). Shown are [CI]/13CO, [CI]/12CO and 12CO/13CO. Color (hne) 
denotes the valne of the line ratio, and the satnration (intensity) its S/N. This velocity 
range corresponds to the Norma far-spiral arm crossing. 
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Fig. A.2.— Evans plots for the velocity range —55 to —40 km/s (left) and their corresponding 
2D color tables (right). Shown are [CI]/13CO, [CI]/12CO and 12CO/13CO. Color (hne) 
denotes the valne of the line ratio, and the satnration (intensity) its S/N. This velocity 
range corresponds to the Scntnm-Crnx near-spiral arm crossing. 

























